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Abstract 

C*"} ■ We present an update of neutral Higgs boson decays into bottom quark pairs 

in the minimal supersymmetric extension of the Standard Model. In particular 
■ the resummation of potentially large higher-order corrections due to the soft SUSY 

breaking parameters and \x is extended. The remaining theoretical uncertainties 
due to unknown higher-order SUSY-QCD corrections are analyzed quantitatively. 
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1 Introduction 



The Higgs mechanism is a cornerstone of the Standard Model (SM) and its supersym- 
metric extensions. The search for Higgs bosons is one of the most important endeavors 
l/" - ) ; at future high-energy experiments. Since the minimal supersymmetric extension of the 

Standard Model (MSSM) requires the introduction of two Higgs doublets in order to pre- 
serve supersymmetry, there are five elementary Higgs particles, two CP-even (h,H), one 
! CP-odd (A) and two charged ones (H^. At lowest order all couplings and masses of the 

Ph| MSSM Higgs sector are fixed by two independent input parameters, which are generally 

chosen as tg/3 = v 2 /vi, the ratio of the two vacuum expectation values Ui )2 , and the pseu- 
doscalar Higgs-boson mass M^. At LO the light scalar Higgs mass has to be smaller 
than the Z-boson mass Mz- Including the one- loop and dominant two-loop corrections 
the upper bound is increased to M h < 135 GeV [T]. The couplings of the various neutral 
Higgs bosons to fermions and gauge bosons depend on the angles a and (3. Normalized 
to the SM Higgs couplings, they are listed in Table [T] 

The pseudoscalar particle A does not couple to gauge bosons at tree level, and its 
couplings to down (up)-type fermions are (inversely) proportional to tg/3. The negative 
direct searches for the Higgsstrahlung processes e + e _ — > Zh, ZH and the associated 
production e + e~ —> Ah, AH yield lower bounds of M h>H > 91.0 GeV and M A > 91.9 
GeV. The range 0.5 < tg/3 < 2.4 in the MSSM is excluded by the Higgs searches for a 
SUSY scale M SUSY = 1 TeV at the LEP2 experiments 1 . 

The scalar superpartners f^R of the left- and right-handed fermion components mix 
with each other. The mass eigenstates /i^ of the sfermions / are related to the current 
eigenstates Jl,r by mixing angles Of, 

fi = hcosOf + f R smO f 

h = ~h sin Of + f R cos Of , (1) 



* This work has been supported in part by the Swiss Bundesamt fur Bildung und Wissenschaft and by 
the European Union under contract HPRN-CT-2000-00149. 

1 The excluded range of tg/3 values depends significantly on the value of the top-quark mass 
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Table 1: Higgs couplings in the MSSM to fermions and gauge bosons [V = W, Z] relative 
to SM couplings. 



which are proportional to the masses of the ordinary fermions. Thus mixing effects are 
only important for the third-generation sfermions t, b, f , the mass matrix of which is given 
by jU 2 
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with the parameters r& = r T = l/r t = tg/?. The parameters Af denote the trilinear scalar 
coupling of the soft supersymmetry breaking part of the Lagrangian. Consequently the 
mixing angles acquire the form 
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and the masses of the squark mass eigenstates are given by 
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The neutral Higgs couplings to sfermions read as jS] 
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(5) 



with the couplings gf listed in Table El 

In this paper we investigate the theoretical status of SUSY-QCD corrections to neu- 
tral Higgs decays into bottom quark pairs. In particular we concentrate on the theoretical 
uncertainties of the partial width in regions, where the SUSY-QCD corrections are large, 
i.e. for large values of tg/? and sizeable magnitudes of the Higgsino mass parameter /x 
jU] . These regions are particularly interesting, since the contributions generated by gluino 
exchange are enhanced by tg/?. They play an important role in the phenomenology of 



2 For simplicity, the £>-terms have been absorbed in the sfermion mass parameters M~ 
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Table 2: Coefficients of the neutral MSSM Higgs couplings to sfermion pairs. 



SUSY-Higgs bosons at high-energy colliders, since they shift the Higgs-boson discovery 
and exclusion regions significantly 0. The corrections can also provide at distinction 
between supersymmetric and non-supersymmetric Higgs bosons. The dominant contri- 
butions have been resummed before [8 . However, the trilinear coupling A^ may be large, 
too. We extend the resummation by including the dominant Aj, terms. 

Although we investigate only the SUSY-QCD corrections, it should be noted that the 
electroweak corrections can be important, too, and yield an additional contribution to the 
uncertainties. The full one-loop electroweak corrections were computed in Ref. 0, and 
later refined in 1 1 GJ including the two-loop contributions to the Higgs boson propagator 
matrix. Section 2 summarizes the present theoretical status of Higgs decays into bottom 
quark pairs and sets the basis for the resummation, which is described in Section 3. In 
Section 4 we analyze the remaining theoretical uncertainties originating from the SUSY- 
QCD corrections in detail for representative MSSM scenarios. In Section 5 we conclude. 



2 Higgs decays into bottom quark pairs 
2.1 QCD corrections 

The partial decay widths of the neutral Higgs bosons $ = h,H,A into bottom quark 
pairs, including QCD corrections, can be cast into the form 

r[* - bb] = ^^ml{M^){gtf [A QCD + Af] . (6) 

where regular quark mass effects are neglected. The large logarithmic part of the QCD 
corrections has been absorbed in the running MS bottom quark mass m^M^,) at the 
scale of the corresponding Higgs mass M$. The QCD corrections A QCD and the top 
quark induced contributions Af read as [TT] 

A QCD = 1 + 5.67^1 + (35.94 - 1.36^) (^MAX 
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where Np = 5 active flavours are taken into account. In the intermediate and large Higgs 
mass regimes the QCD corrections reduce the bb decay widths by about 50% due to the 
large logarithmic contributions. 



2.2 SUSY-QCD corrections 

In the MSSM the full SUSY-QCD corrections to the fermionic decay modes have been 
computed at NLO 0H2I. in the low-energy limit M^, M z ,m b <C m^mj the results can 



be cast into the simple form 



r( 



-> bb) 

K A 

I (a, b, c) 



r 



QCD 



(0 -> bb) 



1 + CpCif,— 

71 

-fy rrig fi tg/3 ^ml^rn^mj) 
1 

1 + 
1 - 



tga tg/3 
tga 



1 + 



tg/3 
1 



ab log — + be log - + ca log - 

6 c a 



(8) 



(a — b)(b — c)(c — a) 

Fqcd(4> — » bb) denotes the QCD-corrected decay width of Eq. (JHJ). It should be noted 
that NLO-terms involving the trilinear mixing parameter Aj, are absent in Eq. (|5|). 



3 Effective Lagrangian and resummation 



3.1 Construction of the effective Lagrangian 

The result of Eq. (jHJ) can be derived from the effective low-energy Lagrangian jS] 1 
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3 This effective Lagrangian has been obtained by integrating out the heavy SUSY particles b, g and 
is thus not restricted to large values of tg/3 only. It should be noted that the scale dependence of the 
running bottom mass and Yukawa coupling is purely QCD-initiated, since the heavy SUSY particles are 
integrated out at a fixed scale of 0{M$usy) and thus do not appear as active partons in the corresponding 
renormalization group equations. 
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after expansion up to NLO. The symbol 4>i{4>2) denotes the neutral components of the 
Higgs doublet coupling to down-(up-)type quarks. The parameter tg/3 = vi/v\ is defined 
as the ratio of the two vacuum expectation values, and v 2 = v \ + v\ = 1/ \f2Gp is related 
to the Fermi constant Gp. The would-be Goldstone field G° is absorbed by the Z boson 
and generates its longitudinal component. The SUSY-QCD corrections turn out to be 
significant for large values of tg/3 and moderate or large /i values. In order to improve 
the perturbative result all terms of O [(a s /itg/?) n ] have been resummed [Hj. The correctly 
resummed effective Lagrangian is given by Eq. 0. The correction Am b is non-decoupling 
in the sense that scaling all SUSY parameters Tn bi2 , nig, /i in Eq. (jlUj) leaves Am b invariant. 
However, its contribution develops decoupling properties P3J, as we will discuss later on. 

Apart from the correction Am b there is a second class of potentially large (non- 
decoupling) contributions at higher orders which may spoil the perturbative reliability 
of the results: The trilinear mixing parameter A b can be of similar size as //tg/3 as e.g. in 
no-mixing scenarios of the sbottom particles. In the low-energy limit of Eq. (jHJ) such terms 
are absent. However, they arise at higher orders. In the following we develop an approach 
to include A b terms in the resummation of Eq. ©. For this purpose we start from the 
unrenormalized effective Lagrangian in the low-energy limit at leading order 



£ L e?f = -W%m + h.C. 



Including higher-order corrections in the low-energy limit, the pole mass m b of the bottom 
quark is given by 



A 

m b = -j=v 1 + T, b (m b ) 

where the self-energy £ D (mb) can be decomposed as 

A A 
S fe (m b ) = [A^i + A 2 v 2 ] = -±vi [Ax + A 2 tg/3] 
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The leading parts in A b and fi are finite at NLO 
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Inserting these two expressions in Eq. (J12)) leads to the well-known result that the radia- 
tive corrections to the bottom mass are proportional to A b — /itg/3, i.e. the off-diagonal 
components of the sbottom mass matrix of Eq. @ . 

The structure of the self-energy beyond NLO can be derived from general arguments 
based on the asymptotic behaviour of the corresponding Feynman- diagrams in the low- 
energy limit. The terms involving A b or \i are generated by mass-insertions in the virtual 
sbottom propagators. At NLO the diagrams of Fig.^ behave asymptotically as 4 
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Figure 1: One-loop contribution to the quantities (a) A x and (b) A 2 . 



a s \ b (A b vi - nv 2 )m g x C o (0, 0; m bi , m^ 2 , rrig) ~ a s m b m g — ' (15) 

M SUSY 

(for M S usy ~ m b ~ Tnt ~ Trig) coinciding with the explicit results of Eq. ()14|). At NNLO 
the leading contributions involving A b and \x are generated by e.g. the diagrams of Fig.El 
The diagrams (a) and (b) behave asymptotically as 



a 2 s \ h {A b vi - /j,v 2 )m g xA (m- bi )x D (0, 0, Ojm^m^mg ,m§) ~ Q? s m h m~ a . ' ^ tg ^ (16) 



while diagram (c) develops the low-energy behaviour 
a%\ b (A b Vi - i2v 2 )mgXB (0]m ii ,mi 2 )xCo(0,0]m' b _,m- b .,mg) ~ a 2 s m b m g Al 2 ^ tg ^ (17) 

Thus, the diagrams of Fig. contribute to the same order as the pure QCD corrections 
to the NLO results and do not generate leading terms of 0(A b ), 0(/i 2 tg 2 /?) nor 0(A b [itg/3). 
This power-counting argument can be applied to all other two-loop diagrams involving fi 
and A b , too. Any further mass-insertion is suppressed by another power of m b /MsusY, 
and is therefore non-leading. 

These arguments can be extended to any perturbative order. Due to the KLN the- 
orem irreducible diagrams do not develop power-like divergences in the bottom 

4 The functions Ao, Bq, Cq, Dq denote the usual one-loop scalar integrals for one-, two-, three- and 
four-point functions. 



6 



1 X 

9 



J — x — ■ — » 

9 



(a) 



(b) 




Figure 2: Non-decoupling two-loop contributions to Ai and A 2 



mass for m b — > 0. Any mass-insertion in the sbottom propagators leads to the replacement 



Therefore, the low-energy behaviour of the mass-inserted diagram is modified by an ad- 
ditional power of m b (A h — fitg/3)/Mg USY . Consequently, the diagrams of Fig. El constitute 
the leading contributions in A b and ntgfl at NNLO. These arguments prove that the re- 
sults of Eq. include all leading powers of a s Ab and a s fitg/3. This is confirmed by the 
explicit two- loop results of Ref . jTHj . 

In order to obtain the effective low-energy Lagrangian from the expression Eq. (jl2j) for 
the bottom mass, we have to perform the replacements v\ — > \f2<p\ and V2 — > V^<p2* i n the 
corresponding bottom mass operator. These replacements lead to the exact interactions 
with non-propagating Higgs fields, i.e. in the low-energy limit of small Higgs momentum 
|17j . The final expression of the effective Lagrangian can be cast into the form 



which differs from previous results by the new factor (1+ Ai) in front of <fi\. This expression 
has to be matched with the renormalized low-energy Lagrangian 
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yielding the relations 5 
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Thus all terms of 0[(a s / Msusy) 
replacement 



1 (fitgf3) m A]!; m ] are resummed by means of the simple 

in the effective Lagrangian of Eq. This proof confirms and extends the resummation 
presented in Ref. jH] and explains the absence of any A b terms in Eq. (|SJ) in terms of 
a clear physical interpretation: the leading A b terms are absorbed in the definition of 
the effective Yukawa coupling X b in the low-energy effective Lagrangian. In a Feynman 
diagrammatic approach this corresponds to a cancellation of the A b terms in the bottom- 
mass counterterms and the genuine irreducible three-point diagrams. This cancellation 
is exact at zero-momentum transfer, but a mild dependence on A b appears when keeping 
all external momenta on-shell due to the momentum dependence of the one-particle- 
irreducible (1PI) three-point functions. 

The final results for the resummed partial decay widths can be cast into the form [see 
Eqs. (PS] 6 
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with the resummed couplings [see Eqs. (l9JlfJJ2()|) ] 
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(23) 



3.2 Validity of the low-energy approximation 

The expression in Eq. (ITT?|) resums the terms of 0[(a s /M S usY) n (^sP) m A^ m ] to all orders 
in perturbation theory. However, there are other kinds of non-decoupling terms in the 
1PI self-energies, as can be inferred already from the NNLO expressions of Eqs. (|16ll7j) . 
The question about the numerical size of these non-leading terms arises, and whether the 

5 It should be noted that the bottom wave-function renormalization constants do not contain any 
leading non-decoupling contribution in Ab and /i. Moreover, it should be emphasized that the combination 
Ai, — /itg/3 only appears in the definition of the bottom mass, while Ab and /itg/3 contribute in a different 
way to the bottom Yukawa coupling and Higgs decay processes. 

6 In order to avoid an artificial singularity in T(h — > bb) for vanishing a the remainder proportional to 
(C - Ci E ) is multiplied by the unresummed Yukawa coupling gf . 
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NNLO resummation is necessary in practical applications. Eqs. ()15H17|) imply that the 

(2) (2) (2) 

irreducible NNLO corrections A\ and A 2 to the selfenergy are of the order of ^{12} ~ 
a s A{ 12 }, while the reducible diagrams contribute as (A{ 12 }) 2 . For the irreducible dia- 
grams to be dominant compared to the reducible ones, the condition (A{ 12 }) 2 < I ^11,2} I ~ 
a s |A{i i2 }| has to be fulfilled, i.e. |A{i j2 }| < a s ~ 0(10%). Therefore, the scenarios with 
the NNLO 1PI being dominant lead to |A|j 2 -.| < (9(1%), so that the NLO corrections are 
small, and the size of the NNLO corrections is of the same order as the deviation of the 
full results from the zero-momentum approximation. This argument can be extended to 
higher orders in perturbation theory. At the n-loop level the non-decoupling 1PI diagrams 
originate from a single vacuum insertion (analogous to the diagrams of Fig. |2J) which are of 
0(a™m\,rng(Ab — fitgP)/Mg USY ) — a" _1 A{i j2 }. Hence, they are negligible, because either 
they are much smaller than the n-loop reducible contribution, or the numerical value of 
the leading corrections is small already at NLO. 

The trilinear mixing parameter Ab cannot be much larger than Msusy, since otherwise 
the color and charge symmetries would be broken ^H]- Thus, the contribution Ai of 
Eq. (j!4|) reaches maximal values of 0(10%), while the term Am;, can be larger by an 
order of magnitude. 

In Fig. El we compare the relative NLO corrections including the resummation of Am.;, 
with the novel NNLO contributions A x of Eq. (jl(Jj) as a function of the pseudoscalar Higgs 
mass Ma for all three neutral Higgs states in the following MSSM scenario with large A^. 

30 

2 TeV 
1.6 TeV 

[A cot (3 
-//tg/3 

-150 GeV (24) 

The relative corrections are normalized to the QCD-corrected decay widths Tqcd(4> ~~ *■ bb) 
of Eq. © in both cases. While the Arrib effects are of 0(10%) and thus of moderate size, 
the novel Ai contributions turn out to be of (9(1%) apart from the small heavy scalar 
Higgs mass range, where they can reach a similar magnitude as the Am;, terms. This 
particular scenario, however, has to be considered as an extreme case. In general the Ai 
terms are small, confirming the previous qualitative discussion. 
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4 Numerical Results 

The numerical analysis of the neutral Higgs boson decays into bottom quark pairs is 
performed for the 'small a e //' MSSM scenario JH] as a representative case: 

tg/3 = 30 

Mq = 800 GeV 

M § = 500 GeV 
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Figure 3: Relative corrections due to (a) the SUSY-QCD corrections including the re- 
summation of Ami, of Eq. Ftty) and (b) due to Ai of Eq. fcTfy as a function of the pseu- 
doscalar mass Ma for all neutral Higgs bosons. The relative corrections are normalized 
to the QCD-corrected decay widths Yqcd{4> ~^ bb) of Eq. (GJ) in both cases. 
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We use the RG-improved two-loop expressions of Ref . |2D] • The bottom quark pole mass 
has been chosen to be Mb = 4.62 GeV, which corresponds to a MS mass TO&(m&) = 
4.28 GeV. The strong coupling constant has been normalized to a s [Mz) = 0.119. 

The resummation effects discussed in the previous section have been derived in the 
low-energy limit M|, M\,m\ <C Mg USY . The question arises, how reliable this ap- 
proximation works in phenomenological applications. In particular, the magnitude of 
0(M^/Mg USY , Mz/Mg USY ,ml/Ms USY ) terms matters for sizeable masses of the low- 
energy particles. This can be tested explicitly by comparing the approximate results of 
Eq. (jHJ) with the full one-loop result. A typical example is depicted in Fig.0]for the 'small 
0£ e ff' scenario, where the relative difference between the full and approximate one- loop 
contributions [see Eq. (jHJ)] 



is presented for all neutral Higgs particles as a function of the pseudoscalar Higgs mass 
Ma- It is clearly visible that the approximation turns out to be sufficient for the heavy 
neutral Higgs particles H, A, but fails for the light scalar Higgs boson h in the decoupling 
limit [21]. However, in the decoupling limit the size of the approximate SUSY-QCD 
corrections strongly decreases, since tga — > — l/tg/3 and thus 



so that the SUSY-QCD corrections become negligible. Due to this behaviour the low- 
energy approximation is sufficient for most phenomenological applications. This also 
explains the failure of the approximation in this case: the large non-decoupling contribu- 
tions from Ami, cancel to a large extent in the lightest Higgs boson couplings, leaving a 
small remainder of the same order as the non-leading contributions. On the other hand, 
this cancellation does not occur for the heavy Higgs bosons, and the effective Lagrangian 
approach yields a good approximation. 

There are two basic sources of systematic uncertainties originating from the SUSY- 
QCD contributions: (i) The MSSM masses and couplings involved in the NLO SUSY- 
QCD corrections will only be known with a sizeable uncertainty at the LHC, while future 
e + e _ linear colliders in the 500 GeV to 1 TeV range will enable precision measurements 
of the SUSY masses and couplings. These errors in the input parameters generate sys- 
tematic uncertainties for the prediction of the partial decay widths, (ii) Due to missing 
higher order results the scale dependence of the strong coupling constant a s will not be 
compensated. The scale variation yields an estimate of the purely theoretical SUSY-QCD 
uncertainty, which will be analyzed quantitatively in this section. 7 

7 The electroweak contributions introduce additional uncertainties, which are not taken into account. 
They provide contributions to Ai and A2 in addition to the SUSY-QCD part. 
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Figure 4: Relative deviations 5^ of the approximate low-energy one-loop result from the 
full NLO expression as a function of the pseudoscalar mass Ma in the 'small a e ff ' scenario 
(a) for all neutral Higgs bosons. For the heavy scalar and pseudoscalar Higgs bosons the 
deviations are negative. The values shown have to be changed in sign. 



The central scale /i of the strong coupling constant appearing in the SUSY-QCD 
corrections will be chosen as the average mass of the involved SUSY particles, i.e. 



m ' bl + m b 2 + 



/<n = ^ ' ' 3 b2 ' " 9 (28) 

In order to estimate the residual scale dependence the scale of a s will be varied between 
/Xo/3 and 3/i - The usual QCD corrections have been included up to the three-loop order 
so that the residual purely QCD-induced scale dependence ranges below the per-mille 
level and can thus safely be neglected. 

The results for the partial decay widths are shown in Fig.Et for the light scalar Higgs 
boson, in Fig. Eh for the heavy scalar Higgs boson and in Fig.Efc for the pseudoscalar Higgs 
boson. These results include the QCD corrections up to NNNLO of Eq. (jHJ) and the full 
NLO SUSY-QCD corrections of Eq. (jHJ) with the resummation of the leading Am& and 
Ai terms according to Eqs. (|9I21|) . It can clearly be inferred from these figures that the 
remaining uncertainties due to the scale choice are typically of the order of 10%. However, 
they are significantly enhanced in regions where the SUSY-QCD corrections become large, 
as in the 'small a e ff '-scenario, which develops a strongly suppressed partial decay width 
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(b) r(H^bB)[GeV] 
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Figure 5: Partial decay widths T(0 — > bb) of (a) the light scalar, (b) the heavy scalar and 
(c) the pseudoscalar Higgs boson in the 'small a e ff' scenario. The shaded bands reflect 
the uncertainties due to the scale choice of the strong coupling constant a s . 
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T(h — > bb) for pseudoscalar masses M4 ~ 150 GeV 8 . This, however, corresponds only to a 
tiny region in the light scalar Higgs mass close to its upper limit for large Ma within 
the 'small a e ff '-scenario. The theoretical uncertainties turn out to be large at M4 ~ 150 
GeV. 

The uncertainties in the partial decay widths T((j) — > bb) translate into systematic 
errors in the corresponding branching ratios. They are depicted in Figs-E^-c for the three 
neutral Higgs bosons. These results have been obtained with the program HDECAY |22] 
after including the results obtained in this analysis. Since the partial decay into bb pairs 
is dominant in nearly the entire Higgs mass ranges, its uncertainty due to the scale choice 
above reduces to a level of 0(1%). However, the scale dependence of T(0 — > bb) develops 
significant systematic errors in the non- leading branching ratios into t + t~, gluon and tt 
pairs. These can reach a level of 0(10%) and are larger than the expected experimental 
accuracy at future e + e~ linear colliders, which clearly calls for a NNLO calculation of 
the SUSY-QCD part. These theoretical errors have to be added to the uncertainties due 
to inaccuracies of the input parameters as presented in j2H] and the theoretical errors of 
the Higgs masses and couplings 3J 9 . They constitute a significant source of uncertainty. 
An analogous analysis is required for the theoretical uncertainties due to the SUSY- 
electroweak corrections beyond NLO. However, this is beyond the scope of our paper. 



5 Conclusions 

In this paper we have reanalyzed the neutral scalar Higgs decays into bb pairs in the MSSM 
with particular emphasis on the SUSY-QCD corrections and their theoretical uncertain- 
ties. We have extended the resummation of large non-decoupling SUSY-QCD corrections 
of C(a s /itg/3 / 'M S usy) by the inclusion of non-decoupling terms of 0(a s A b /M SUSY ) which 
have not been taken into account in previous analyses. We have shown that these terms 
are absent at NLO in the effective Lagrangian but arise at NNLO and beyond. This can 
easily be traced back to the renormalization of the bottom Yukawa coupling in the low- 
energy limit, where the heavy SUSY particles are integrated out. We have obtained the 
important result that these novel contributions hardly affect the theoretical predictions 
for the partial decay widths into bb pairs so that they do not endanger the reliability of 
the perturbative result in contrast to the leading terms of 0(a s fj,tgP / Msusy) ■ 

We investigated the remaining theoretical uncertainties generated by the SUSY-QCD 
corrections quantitatively. While the theoretical errors of the partial decay widths — > 
bb) turn out to be of 0(10%), this effect cancels to a large extent in the branching ratios 
BR(<p — > bb) due to its dominance. It appears, however, as a sizeable increase in the sys- 
tematic uncertainties of the non- leading branching ratios into r + r", gluon and tt pairs, 
which appear to be larger than the anticipated experimental accuracies at future linear 
e + e~ colliders. This clearly calls for a NNLO calculation of the SUSY-QCD part, which 

8 The explicit value of the pseudoscalar mass where the Yukawa coupling vanishes depends strongly 
on the included higher-order corrections. 

9 The uncertainties due to the Higgs masses will be eliminated to a large extent, once they will be 
measured directly in future experiments. 
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Figure 6: Branching ratios of (a) the light scalar, (b) the heavy scalar and (c) the pseu- 
doscalar Higgs boson in the 'small a e ff ' scenario. The shaded bands reflect the uncertain- 
ties due to the scale choice of the strong coupling constant a s . 
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is beyond the scope of this work. 
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